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Sliding Discrete Fourier Transform (SDFT) is very effieient regarding 
computational load and it possesses a very fast phase angle detection with 
excellent harmonic rejection at nominal frequency. However, at off-nominal 
frequency, SDFT generates errors in both magnitude and phase angle due to 
spectral leakage. This paper introduces a workaround for Fourier Transform 
to handle this disability under off-nominal frequency while avoiding 
variable-rate sampling. Sliding Fourier Transform (SFT) is used as a phase 
detector for a phase-locked loop whose output frequency is used to drive the 
SFT. The paper revisits the mathematics of Fourier Transform (FT) in a 
three-phase setting via a time-domain approach to show a newly proposed 
filtering technique for the double-frequency oscillation just by summing the 
FT sine/cosine filter outputs of the three individual phases. Also, the analysis 
aims to determine and correct the phase and magnitude errors under off- 
nominal frequency operation. The proposed technique (SFT-PLL) is tested in 
real time on dSPACE DS1202 DSP using voltage vectors that are pre¬ 
generated to simulate the most adverse grid conditions. The testing scenarios 
compare the performance of the SFT-PLL with that of the Decoupled 
Stationary Reference Frame PEL (dapPLL). The results prove that SFT-PLL 
is superior to dapPLL. 
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1. INTRODUCTION 

Fuel eells and other Renewable Energy Sourees (RES), sueh as PV systems and Wind Energy 
Converters (WECS), use inverters to injeet the eaptured energy into the utility grid [1]. The sueeessful 
integration of these sourees to the eleetrieal grid is still having a lot of issues and teehnieal ehallenges [1]. 
One of the most eritieal issues of these integrations is the synehronization of the grid-side eonverter (GSC) to 
the eleetrieal grid [2]. The synehronization of the GSC to the eleetrieal grid ean be defined as the proeess of 
minimizing the differenee in phase angle, frequeney and voltage between the GSC output and the utility grid 
waveform [3]. This synehronization proeess has to be realized before the eonneetion between the RES or the 
GSC and the utility grid takes plaee, then it is to be maintained after the eonneetion is established. 

In [4], the ideal synehronization teehnique must instantly respond to any grid ehanges, effieiently 
traek the phase angle of the eleetrieal grid, eompetently deteet any variation in the grid frequeney and 
eompetently eliminate harmonie eomponents and disturbanee from the fundamental waveform. There are lots 
of methods used for traeking the phase angle of the grid voltage sueh as zero-erossing deteetion method 
(ZCD), Spaee Veetor (SV) method, Kalman filter, phase-loeked loop methods (PEE) and Diserete Fourier 
Transform (DFT). The problem is none of these methods eompletely satisfies the full requirements of an 
ideal synehronization teehnique 
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Zero-Crossing Detection (ZCD) is the simplest angle synchronization method. The zero-crossing 
point is detected twice every single cycle and the interval between each two successive crossings is measured 
to any desired resolution. This method is very simple but it can be used only if the input signal is a stable 
sinusoidal signal as it is very sensitive to harmonics and noise which can cause multiple contiguous 
zero crossings [5]. 

Phase locked loop is a well-established technique and there are a lot of implementations of PEL 
synchronization. Each implementation differs from the other from the point of view of simplicity, immunity 
to disturbance, robustness in case of grid faults, fast convergence to the grid angle, and how it deals with 
disturbances. Synchronous Frame PEL (SF-PLE) technique is widely used in three-phase grid-connected 
power converters for its simple implementation, fast and accurate estimation of the phase angle of the grid 
[6]. Under ideal grid conditions, without any disturbances or harmonics, this technique exhibits fast-tracking. 
However, it is not a good solution in the presence of unbalanced grid voltages [7]. 

Real and imaginary power phase-locked loop (PQ-PLL) is a robust digital PEE algorithm based on 
the instantaneous real and imaginary power. This algorithm is introduced to maintain synchronization in 
presence of harmonics and negative sequence unbalances. Another algorithm which can deal with unbalanced 
grid voltages is the Double Synchronous Frame PEL (DSF-PLL). It is a very suitable solution for grid- 
connected converters that operate in unstable frequency grid and unbalanced conditions [8]. 

Sinusoidal signal integrator PEL algorithm (SSI-PEL) has the advantage of operating well under 
unbalanced and distorted grid [8]. Another advantage of this algorithm is that it can be used for the single¬ 
phase system with a small modification. Enhanced phase-locked loop (EPEE) is a frequency-adaptive 
nonlinear technique. EPEE shows a high degree of robustness to the harmonic, noise, and unbalanced grid 

[9] . This algorithm can be used for single phase applications. Quadrature PEL (QPLL) is based on estimating 
the in-phase and the quadrature phase value of the fundamental component of the input voltage signal. This 
method is a good solution for distributed generation [10]. Predictive Phase-Locked Loop (PPLL) is a 
complex synchronization method. However, PPLL has a wide-range synchronization capability and the 
synchronization information can be taken out within two voltage cycles under the worst case. Stationary 
Frame PEL (a|3PLL) in [11], Decoupled Double Synchronous Reference Frame PEL (ddsrfPLL) and 
Decoupled Stationary Reference Frame (da|3PLL) techniques have the widest use in industry. The limitation 
of the apPLL is its disability to track the grid angle in the case of unbalanced faults. The ddsrfPLL technique 
overcomes this limitation by decoupling the negative and the positive sequence of the grid voltage but the 
severe drawback of this technique is the high overshoot in the estimated grid angle when the fault occurs. 
The dapPLL overcomes this drawback but it spends more than two cycles to track the grid frequency in case 
of frequency change. 

Discrete Fourier transform method (DFT) is one of the earliest methods used for phase angle 
estimation. This method provides a high degree of immunity against harmonics [12]. Sliding Discrete Fourier 
Transform (SDFT) is an improved version of the DFT in which optimized use of computation power and 
reduced latency are provided by using a recursive formulation of the DFT algorithm. However, if the time 
window of the SDFT is not identical with the grid period, a phase shift occurs between the estimated angle 
and the actual grid angle. In [12] some strategies have been reported to compensate for this phase error such 
as the use of time-varying Fourier coefficients and polynomial regression analysis. In fact, these strategies 
are rendered with noticeable complexity. Other fixing methods for the SDFT drawback are presented in [13]. 
However, the presented methods involve the use of a PI controller and a Numerically Controlled Oscillator 
(NCO). The method adopts a variable sampling rate of the grid voltage controlled by the NCO. Variable rate 
sampling adds difficulty to embedding the SDFT algorithm in the main application DSP controller which 
normally runs at fixed execution rate. Although the two proposed methods in [13] are simpler than those of 

[10] , both methods compensate the phase error on a cycle-by-cycle basis, which is rather a slow process. The 
undesirable features of the SDFT in the case of off-nominal frequency are not only limited to the phase error 
which is a function of the frequency deviation but also includes a high-frequency oscillation observable in the 
magnitude and phase angle estimates. This high-frequency oscillation is referred to in the literature as 
double-frequency ripple and some post-processing for frequency and phase angle estimates in the off- 
nominal operation mode are proposed. The post-processing can be done either by simple digital averaging 
filter or by resampling filter which is relatively more complex in implementation. 

From the previous survey, it appears that there is still some room for improvement in the phase 
angle tracking issue especially in the case of variable and off-nominal grid frequency accompanied by high 
harmonic content and/or in unbalanced three-phase systems. This paper revisits the FT concept through a 
time-domain analysis to further analyze the undesirable features associated with the off-nominal grid 
frequency. An efficient and natural filtering method for the double-frequency oscillations is devised. The 
laws relating the phase and magnitude errors to the frequency difference are derived. The formulation of the 
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FT in a sliding mode is explained. Experimental testing is used to eonfirm the validity of the 
proposed teehniques. 

This paper is organized as follows: the mathematieal analysis for determining the phase and 
magnitude errors under off-nominal frequeney eondition is presented in seetion 2. The reeursive 
implementation is explained in seetion 3. The strueture of the proposed SFT-PLL is explained in seetion 4. 
The method used for experimental testing and the test results in real-time on a DSP prototyping platform 
(DS1202 by dSPACE™) are presented in seetions 4 and 5, respeetively. Seetion 6 presents the eonelusions of 
this paper. 


2. MATHEMATICAL ANALYSIS 

The grid voltages are represented by a balaneed set of three-phase sinusoidal voltage Vg with a 
frequeney fg (Hz) and with a peak value E (Volt), thus Vg is given by: 

Vga=E sm{2nfgt + do), 

Vgf,=Esm{2nfgt-'^ + eo') (1) 

Vgc = E sin(2nfgt -h 2n/3 + 6 q) 

where 6q is the initial phase shift (rad) of the voltage Vg^, this phase shift depends on the random instant at 
whieh the measurement is started (^=0). For the sake of this analysis fg will be assumed eonstant. 

The Fourier transform is based on applying two orthogonal filters, namely a sine filter and a eosine 
filter on the input signal. The point here is that we set the phase sequenee and the phase shift between the 
sine/eosine filters in a way that serve additional filtering purpose as will be seen later on. Therefore, the three 
grid voltages will be proeessed using two sets of balaneed three-phase signals V^x and V^y having the same 
phase sequenee as the grid voltage, running at a pre-known nominal grid frequeney (Hz) and having a 
unity peak value. For the sake of this analysis f^ will be assumed eonstant [14]. 

The first set of signals V^x is given by sine funetions and is ealled the sine or direet set: 

Vcax = sin(27r/„t), 

Vcbx = sin (2nfnt - (2) 

Vccx = sin(27r/„t + 2n/3) 

The second set of voltages V^y is given by cosine functions and is called the cosine or quadratic set: 

Vcay = cos(2nf„t), 

Vcby = cos ^27r/„t - y) (3) 

Vccy = cos(2nfjit -h 2n/3) 


The Fourier transform requires performing the following integrations to get the angle 9q\ 


L-r I'f 


^cix- 



^ciy 


( 4 ) 

( 5 ) 


(4) and (5) are eompaet forms for the three-phase where the subseript i = a, b, c and 7^ = 1/A- For 
any positive and non-zero value of A? two eases are eonsidered for ealeulating these integrals, the first ease is 
fg^fn (off-nominal frequeney ease) whieh is the generie ease and the seeond ease fg=fn (nominal 
frequeney ease) whieh is a speeial ease of the first ease. 
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2.1. The first case: fg ^ /„ (off-nominal frequency case) 

¥ or \fg — f^\ < , and using suitable trigonometric identities it can be proven that in (4) (let us 

call it the direct component) and in (5) (let us call it the quadratic component) are functions of time and 
frequencies fg , f^ that are given as follows: 

For phase a: 

^a(0 = COs[a)£)t + ^0 + “ /C2COS[u)7’t + ^0 + ^^ 2 ] (6) 

ya(0 = sin[a)£,t + ^0 + ^^ 1 ] + k 2 sm[a)Tt + 9 q + A 62 ] (7) 

For phase b: 

x^(t) = /ci cos[a)j)t + 6 q + AOi] — k2Cos[a)Tt + 6 q + A 62 + 2pi/3] (8) 

y^(t) = /ci sin[u)£,t + + ^^ 1 ] + /c2sin[u)7^t + ^0 + ^^2 + 2pi/3] (9) 

And for phase c: 

Xcit) = /ci cos[a)£,t + ^0 + ^^ 1 ] “ /c2Cos[a)7’t + + ^^2 “ 2pil3\ (10) 

yc(0 = sin[a)£,t + ^0 + ^^ 1 ] + k2sm[a)Tt + 9 q + A92 - 2pil3] (11) 

where 


Wd = 2n{fg - fn) 

( 12 ) 

"r = 27r(4 + fn) 

(13) 

> 

II 

(14) 

^ E sin(A0i) 

“ 2 A6>i 

(15) 

OJj' 

Ado = — 

" 2fn 

(16) 

^ E sin(A 6 l 2 ) 

' “ 2 A 0 , 

(17) 




Figure 1. Waveforms (9/Xi(t),Xt(t),yi(t),yt(t): F = lV,fg = 60Hz, ^ = 50 Hz and “'n:/2 
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It can be noticed that both of the direet and quadratie eomponent Xj and yi for the three-phases (i = 
a, b, c) as shown in (6) to (11) are eomposed of two sinusoidal terms. One of them is larger in magnitude and 
its frequeney is the differenee between the grid aetual and nominal (or filter) frequeney (if the two 
frequeneies beeome equal, then it is a DC value as will be shown in (6), (7) of the seeond ease). The other 
has a mueh smaller magnitude and its frequeney ( 0 )^) is the sum of the grid aetual and nominal frequeney (if 
the two frequeneies beeome equal, then it vanishes as the magnitude beeomes zero as will be eonfirmed in 
the seeond ease). Referring to the seeond term in the right hand side of eaeh of (6), (8) and (10), it is evident 
that the three terms eonstitute a balaneed three-phase set in eosine notation whose sum is zero. Thanks to the 
proper seleetion of the sine and eosine filter phase angles. Referring to the seeond term in the right hand side 
of eaeh of (7), (9) and (11), the three terms also eonstitute a balaneed three-phase set in sine notation whose 
sum is zero. Here it beeomes elear the usefulness of the summing proeess eombined with the proper seleetion 
of the FT filter phase angles. The summing proeess thus results in filtering the double-frequeney oseillations 
without any additional or post-filtering. Some representative values (E = l,fg = 60 Hz, = 50 Hz and 
6q = — 7i/2) are used to plot 

Figure 1 where the waveforms of Xj and are plotted in the left subplot and the right subplot 
respeetively against time to illustrate the double-frequeney ripple, the sums x^ and y^ are also plotted in the 
two subplots to show the smooth sinusoidal variations as the three higher frequeney oseillations running at 
ifg + fn) caneel eaeh other. Therefore, the sum of the direet x^ and the quadratie eomponents y^ for the three 
individual phases, are given by (18) and (19), respeetively: 


Xt = ^ Xj = 3 kiCos[a)Dt + Oq + AO^] 

i=a 

c 

Yt = ^ Yi = 3 kiSin[(ODt + 00+ Aej 


Dividing (19) by (18) yields: 

— = tan[cL)Dt + 00 + A0i] 
Therefore: 

27i(fg - fn)t + 00 + A01 = tan“^ 


(18) 


(19) 


( 20 ) 


( 21 ) 


The phase angle of the grid voltage Vg^ that is equal to 2nfgt + 6 q from (1) is then given by 
rearranging (21) as follows: 


0a = 27ifgt + 00 = tan ^ + 27ifnt - A 01 (22) 

Substituting from (14) into (22) yields, 

e,=tan-i(j)+ (23) 

Reealling that fg is eonstant, then ^ = 0 and A6^ is eonstant. 6 q is also eonstant. Therefore, the 
time derivatives of both A^^and ^o are typieally zeros. Consequently, the value of fg ean be ealeulated by 
taking time derivative of both sides of (22) as follows: 

Onee fg is estimated, A^^is determined and eonsequently ^^ean be estimated from (22). To estimate 
the magnitude E, it ean be found from (15), (18) and (19) that: 
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E 


2 A01 

3' sin(A0i)' 




(25) 


sm(A0i) 

called the uncompensated magnitude. 


The term ^ can be called the magnitude correction factor and the term ^ + yt^ can be 


2.2. Second case: fg = (nominal frequency case) 

By using suitable trigonometric identities, it can be proven that Xj and yi are independent of time 
and are given as follows: 


E cosC^o) 

Xa=Xi,=Xc= 2 


(26) 

E sin(0o) 

ya = yb=yc= ^ 


(27) 

Therefore, the summing for x andy components of the three-phase are given by (28) and (29), respectively: 

^ 3 E cos(0o) 

i=a 


(28) 

^ 3E sin(0o) 

yt=2,3^^= 2 

i=a 


(29) 

Now dividing both sides in (29) by their respective sides in (28) yields: 



yt _ sin(0o) ^ 

Xt cos(0o) “ 


(30) 



(31) 

Recalling the equality of ^and fg, the phase angle 6a of grid voltage Vga 

is then given by (32): 


6a = 2nfgt + 6q = 2nfat + tan~^ 


(32) 

According to (23), the last term in the right-hand side of (23) vanishes when fg=fn thus leading to 
(32) and confirms that the nominal frequency case is a special case of the off-nominal frequency case. 
Therefore, it is obvious that (23) can be generalized for both nominal and off-nominal frequency cases. Thus, 
phase angle tracking is a seamless process when the actual frequency shifts from nominal to off-nominal 
mode and vice versa during actual operation. 

The magnitude E can simply obtained from (28) and (29) as follows: 

7 

^ = o-+yt^ 


(33) 


A0 

Knowing that: lim —= 1, the (25) becomes the same as (33) and it turns as stated earlier that 

the nominal frequency case is a special case of the off-nominal generic case. Therefore (25) can be a valid 
expression for both cases. Thus, the magnitude tracking is also a seamless process should the actual grid 
frequency shifts from off-nominal to nominal mode of operation and vice versa. 

Now it is clear that by proper selection of the phase sequence and phase shift of the waveforms used 
for sine and cosine filters of the FT and applying the summing procedure we get a very effective filtering for 
the double-frequency oscillation characterizing the FT in off-nominal frequency case. Figure 2 illustrates the 
schematic diagram of the proposed FT implementation based on (1) through (33). In the leftmost block, the 
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sine and cosine filters are applied to the three-phase voltage, and the summing procedure is applied to get 
and y^. The uncompensated phase angle given by (21) is estimated in the leftmost block as well. The fg is 
estimated according to (24) in the intermediate block. Using the estimated fg, the phase angle correction term 
A6i can be estimated in the rightmost block. Subtracting the correction term from the uncompensated phase 
angle as per (23) results in the true phase angle of the grid voltage Vga- The magnitude correction factor can 
be easily calculated also from (25) in the bottom block and finally the compensated magnitude can 
be estimated. 


In fgt H- A 01 + Oq 



Figure 2. Schematic diagram of the proposed FT implementation for three-phase system 


3. SLIDING FT IMPLEMENTATION (SET) 

The sliding integrations from ^ to t -I- 7^ as mentioned in (4) and (5) can be carried out numerically 
in a discrete sampled system. The method applied for obtaining and updating the numerical integration over 
one cycle upon picking each fresh sample of xi (or yi) is to insert a delaying buffer of length 7^ after the 
discrete time integrator [15]. Thus, integration over 7^ is updated every new sample of xi and it is 
implemented as shown schematically in Figure . In this implementation, the buffer is initialized with zeroes 
and it takes 7^ 77^ samples to be filled incrementally with valid readings as 7^ is equal to the sampling period. 
Based on this, the detected grid angle and frequency needs one complete cycle (7^) to converge to their 
actual values for all ^ ^ f^. It can be shown that for fg = the detection is almost instantaneous. 



t+Tn 


X. dt 


Figure 3. Updating Numerical integration over one 7^ cycle upon each new sample of x 


4. THE PROPOSED SFT-PLL 

The implementation of SFT shown in Figure 2 works very fine and gives fast and accurate results in 
the nominal and off-nominal frequency cases as long as the input voltages are free from harmonic distortion 
Once the input frequency shifts from nominal value while harmonics are existing the estimated phase angle, 
the estimated magnitude and the estimated frequency all suffer from oscillations due to spectral leakage and 
the SFT loses its excellent harmonic rejection and becomes unpractical. This is shown using numerical 
simulation in Figure 4. Figure 4 (a) illustrates that the input voltage is distorted with harmonics starting from 
t=t)A sec to /=0.45 sec and the frequency is stepped up then down by 10 Hz at ^0.15 and /=0.35 sec 
respectively. Figure 4 (b) and (c) show that the SFT rejects the harmonics after exacty one power cycle and 
gives accurate frequency and phase angle estimation when it works at nominal frequency (50 Hz) i.e. from 
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/=0.05 sec to /=0.15 sec and from /=0.35 sec to /=0.5 sec. However, during the off-nominal period from 
/=0.15 sec to /=0.35 sec, the existence of harmonics characterizes the estimation of both frequency and phase 
angle by a sustained error. 



time (sec.) 


(a) 



(b) 



Figure 4. (a) Voltage of phase a (other phases omitted for clarity), (b) estimated frequency by SET (c) error in 

phase angle estimation by SET 


To overcome this deficit, and considering the ability of the FT to detect the phase difference Oo 
between the angle of the sine/cosine filters and the angle of the input voltage as shown by (23), the SET can 
be used as a phase detector in a standard phase locked loop. To do this, the SET time window (or number of 
processed samples) should be continuously and properly adjusted to enforce nominal-frequency operation a n 
and hence maintain its excellent harmonic rejection. This is schematically shown in Figure 5 and Figure 6. 
Note that the shaded rectangles in phase b and c in Figure 5 are similar in construction to that of phase a. The 
phase-locked loop PID controller in Figure 6 sets the loop frequency so as to eliminate the phase difference 
A6. The frequency fg and the phase angle 6 g estimated by the PEE are fed back to the variable window width 
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SFT to calculate the sine and cosine filters whose outputs are used to explore the phase difference A9. The 
PLL therefore adjusts the frequency of the SFT such that when it locks, the SFT works at the actual grid 
frequency. In this case, the filter frequency shifts to fg and the operation mode shifts from off-nominal to 
nominal and A6 vanishes. This arrangement provides fast phase difference detection while maintaining the 
SFT well-known excellent harmonic rejection. Figure 5 and Figure 6 together illustrate the strueture of the 
proposed SFT-PLL teehnique. 



Figure 5. Sehematie representation of the proposed variable window width SFT phase deteetor 


/n 



o 


Figure 6. Phase-loeked loop sehematie representation 


5. EXPERIMENTAL TESTING METHODOLOGY & RESULTS 

The experimental testing insures the applieability of the proposed teehnique in real-time setting and 
under praetieal operating eonditions that may violate some of the eonsidered assumptions made during the 
mathematieal analysis. The testing method depends on generating the voltage veetors that eontains the target 
test eharaeteristies (DC offset, harmonie distortion, magnitude, frequeney or phase jump and unbalanee of the 
three-phases) using offline MATLAB simulation. These voltage veetors are sampled at 64 sample per power 
eyele (0.02 see) and reeorded along with an array of eorresponding timestamps. Also the true phase angle, the 
true frequeney and the true magnitude whieh areused for generating them are sampled and reeorded for sake 
of eomparison and validation. Then these voltage veetors are played baek in real-time at the same sampling 
rate using dSPACE dsl202 DSP D/A ehannels, read baek via A/D ehannels and are used as inputs to the 
SFT-PLL. The output of the SFT-PLL are then reeorded in MAT files, then eompared with the true values 
that are pre-reeorded to eompute the errors. The same proeedure is repeated with da|3PLL teehnique and the 
error in both eases are plotted on same axis to eompare the performanee of both teehniques. The testing 
voltages are: 

a. 50 Hz balaneed three-phase (peak E = 220 V) with the following testing seenarios: Harmonie distortion 
(eomposed of third, fifth, seventh and ninth harmonie eomponents with 70%, 60%, 30% and 20% of the 
fundamental magnitude respeetively) switehed-on at t^l see and switehed- off at ^1.5 see. Phase jump 
of nil at see, frequeney step-up of 5 Hz at t^2.5 see. Magnitude step- up of 0.2 E at ^3.5 see. 
Adding DC offsets to balaneed three-phase at t=A see, 0.1 E to phase A, -.05 E to phase B, -0.05 E to 
phase C. 


Int J Pow Elee & Dri Syst, Vol. 10, No. 2, June 2019 : 1034 - 1048 









































































Int J Pow Elec & Dri Syst 


ISSN: 2088-8694 


□ 1043 


b. Unbalanced three-phase input with knowing symmetrical components structure as follows: magnitude 
of positive sequence is 220 V and its angle is zero, magnitude of negative sequence is 60 V and its angle 
is -150 deg., magnitude of zero sequence is 20 V and its angle is 10 deg. 

In all of these tests, the error is computed as the difference between the true quantity and the 
estimated one. In the case of unbalanced three-phase, the true value for phase angle estimation error is the 
phase angle of phase A of the positive sequence. 

Harmonic distortion test: Figure 7 (a) shows the introduction of harmonics at t=l sec. Figure 7 (b) 
and (c) show that da|3PFF gives continuous oscillation in estimation of both frequency (±20 Hz) and angle 
(±4 deg.) while SFT-PFF gives oscillation in estimated frequency (±lHz) for only one power cycle. 
Oscillation in estimated angle continues for one power cycle and error settles to zero in four cycles. 



(a) 



(b) 



(C) 

Figure 7 (a) Voltage of phase A, harmonics switched on at ^1 sec, (b) frequency estimation error, harmonics 
switched on at /=1 sec, (c) phase angle estimation error, harmonics switched on at /=1 sec 
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Phase jump test: In this test a izH phase jump at t=2 see is done. Figure 8 (a) and (b) show that 
dapPLL frequeney error peak is 28.6 Hz. The frequeney error settles to zero in more than 6 power eyeles. 
Phase angle error settles to zero in more than 6 power eyeles. SFT-PLL results in frequeney error peak of 
13.5 Hz. The frequeney error settles to zero in less than 5 power eyeles. Phase angle error settles to zero 
within 5 power eyeles. Henee, SFT-PLL gives superior performanee. 



tlm« (sac.) 


(a) 



(b) 

Figure 8 (a) frequeney estimation error, tiH phase jump at t^2 see, (b) phase angle estimation error, tiI2 phase 

jump at t=2 see 


Frequency step change test: Figure 9 (a), and (b) show that for da|3PLL, the frequeney estimation 
error settles to zero in more than 8 power eyeles. Peak phase angle error is 8.2 deg. The error in phase angle 
settles to zero in more than 10 power eyeles. For SFT-PLL, the frequeney estimation error settles to zero in 5 
power eyeles. Peak phase angle error is 12.8 deg. Error settles to zero in 5 power eyeles. Therefore, SFT-PLL 
errors settle faster than those of da|3PLL. 



(a) 


Int J Pow Elee & Dri Syst, Vol. 10, No. 2, June 2019 : 1034 - 1048 

































Int J Pow Elec & Dri Syst 


ISSN: 2088-8694 


□ 1045 


Figure 9 (a). Frequency estimation error, 5 Hz frequency step up at t=2.5 sec 



(b) 

Figure 9. (a) frequency estimation error, 5 Hz frequency step up at t=2.5 sec, (b) phase angle estimation error, 

5 Hz frequency step up at t=2.5 sec 


Magnitude step change test: In this test, a 20% step change in magnitude is done at /=3.5 sec. 
Figure 10 (a) and (b) show that da|3PFF results in decayed oscillation in frequency with a peak -0.75 Hz and 
lasting for 4 power cycles. Phase angle error of 1 deg. decays to zero in 4 power cycles. While for SFT-PFF 
there is no error at all in frequency or phase angle. SFT-PFF is thus more robust against magnitude 
disturbances than da|3PFF. 



(a) 



(b) 

Figure 10. (a) Frequency estimation error, 20% step up in magnitude at t=2.5 sec, (b) phase angle estimation 

error, 20% step up in magnitude at ^3.5 sec 


DC-offset test: In this test, the three phase voltages are affected with the prescribed DC offsets at 
/=4 sec. Figure 11 (a) and (b) show that da|3PFF results in continuous oscillation in frequency error of ±2.2 


Frequency adaptive sliding fourier transform for synchronizing VSI to the grid (Osama M. Arafa) 



















1046 □ 


ISSN: 2088-8694 


Hz and in phase angle error of ±2 deg. While for SFT-PLL the oseillation in frequeney error and phase angle 
error has a smaller magnitude and notieeable only within two and four power eyeles respeetively. SFT-PLL 
is thus more robust against DC offset disturbanee than da|3PLL. 



(a) 



4 *.02 4.04 4.BB 4.Da 4.1 4.1 Z 4.14 
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(b) 

Figure 11, (a) Frequeney estimation error, three-phase voltage with DC offset (Va+20, Vb-lO, Ve-lO), (b) 
phase angle estimation error, three-phase voltage with DC offset (Va+20, Vb-lO, Ve-lO) 


Unbalanced three-phase test: Figure 12 (a) shows the three-phase voltages with the preseribed 
unbalanee. Figure 12 (b) and (e) show that da|3PLL results in eontinuous oseillation in frequeney error (-2.9 
Hz:+2.6 Hz) and in phase angle error of ±3 deg. For SFT-PLL, the deeays to zero in 4 power eyeles. While 
for SFT-PLL there no errors at all in estimating frequeney and phase angle. SFT-PLL is thus is eompletely 
immune against unbalanee disturbanee while dafBPLL is not. Note that the referenee phase angle in this ease 
is the positive sequenee phase angle. 
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time (sec.) 


(b) 



(c) 

Figure 12. (a) three-phase unbalanced set, (b) frequency error in case of three-phase unbalanced input (c) 
frequency error in case of three-phase unbalanced input 


6. CONCLUSION 

The mathematical analysis in time domain of the SET of a single frequency three-phase set shows 
that by properly setting the phase angle used in the sine/cosine filters of the individual phases and then by 
summing the direct components and the quadratic components independently, the double-frequency 
oscillation is completely removed without post filters. The phase angle error and the magnitude error are also 
analyzed and their values are determined analytically. Therefore, any errors in the phase angle and magnitude 
can be corrected. However, due to spectral leakage, the SET shows excellent harmonic rejection when 
working at nominal frequency while it shows inferior performance at off-nominal frequency mode of 
operation when the input has a high harmonic content. The benefits of the excellent harmonic rejection of 
SET at nominal frequency mode of operation are preserved by timely adjusting the SET frequency using a 
phase-locked loop. The proposed SFT-PEE has a fixed sampling rate which is consistent with the design of 
common digital controllers. The proposed SFT-PEE has been tested and validated using DSP in real time. A 
comparative study with one of the best performing PEE methods da|3PEE under different practical conditions 
has been presented. The comparison demonstrates that SFT-PEE can successfully handle unbalanced three- 
phase input and computes the positive sequence phase angle while da|3PEE gives a sustained oscillatory error 
in this case. The SFT-PEE is very resistant to magnitude disturbances while da|3PEE is affected by such 
disturbances. Indeed, the SFT-PEE outperforms the da|3PEE in the presence of harmonic distortion or/and 
DC offset. In a very wide range of disturbances, the SFT-PEE error settles to zero faster than that of the 
dapPEE. The excellent performance strongly suggests the SFT-PEE for synchronizing renewable energy 
sources with the utility grid under challenging grid conditions like fault-ride-through periods. 
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